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I .  INTRODUCTION 


A.  OVERVIEW 


and  applied  statistics,  e.g.  in  reliability  and  maintaina¬ 


bility  studies,  is  that  o£  predicting  from  available  data 
the  probability  that  a  future  observation  exceeds  a  given 


value.  An  example  arising  in  nuclear  plant  reliability  is 
that  a  crucial  repair  or  down  time  exceeds  h  (=4.)  hours. 
Another  problem  is  to  predict  the^lOO-year  £lood~^~ox  earth¬ 
quake,  etc.  The  latter  problem  is  difficult  because  there 


will  usually  be  far  less  than  100-years  worth  of  data  to  work 
with.  Still  Mother  problem  is  that  of  predicting  the  pro¬ 
bability  of  survival  for  h  *  5  years  for  a  cancer  victim  re- 

_ .■>-* 

ceiving  a  particular  treatment.  — - - — - -  ~~  " 


The  simplest  formulation  is  to  assume  that  the  data  is  a 

V 

random  sample  from  a  probability  distribution  Fx(x)  (con¬ 
tinuous,  i.e.  having  a  density),  with  density  f  (x) .  That 
is,  observed  values  are  x^,  x2,  .  .  .  ,  xfl,  being  independent 
realizations  of  independent  identically  distributed  random 
variables  generically  denoted  by  X.  If  one  is  willing  to 
assume  also  that  the  mathematical  form  of  F  (x)  *  F  (x;0)  is 
known  (e.g.  is  Log-normal,  or  Gamma,  or  another  candidate) 
then  what  one  can  do  is: 


(a)  Estimate  the  possibly  multidimensional  parameter  9 

2 

(e.g.  9  could  be  y,  a  ,  the  population  mean  and  variance  for 
a  log-normal  model,  estimated  by  lnx  ■'  y  and 

Slnx2  =  1  (lnx^  -  lnx)2  classically). 

i=l 

(b)  Quote  the  point  estimate  Fx(x;9),  or,  in  the  present 
case  1  -  Fx(h;§)  for  probability  of  survival  beyond  h. 

(c)  Utilize  facts  about  the  sampling  distribution  of  § 
to  find  a  standard  error  or  confidence  limits  on  F  (h;9),  the 
survival  probability. 

The  basic  assumption,  then,  is  that  data  can  reasonably 
be  assumed  to  be  a  random  sample  from  a  fixed  distribution, 
the  form  of  which  is  known.  There  are  various  ways  in  which 
such  convenient  assumptions  can  be  violated,  one  obvious  one 
being  that  the  fixed  distribution  idea  is  not  justifiable 
(perhaps  because  of  important  detectable  variation  in  the 
distribution  from  location  to  location,  or  plant  to  plant, 
from  repair  crew  to  repair  crew,  etc.) .  Another  might  be 
that  some  data  points  are  missing:  too-short  ones  (down 
times)  are  not  written  down  or  else  are  recorded  incorrectly, 
and  too-long  ones  are  regarded  as  being  so  exceptional  as 
never  to  recur,  and  hence  are  removed.  Possible  or  likely 
departures  from  the  basic  assumption  should  be  investigated. 
The  raw  data  should  be  carefully  examined  in  an  exploratory 
spirit  (see  J.  W.  Tukey  (Ref.  1]),  e.g.  by  graphics  to  check 


for  departures  from  the  basic  "stationary"  assumption.  In 
this  discussion  we  rule  out  such  variations. 


This  paper  gives  an  account  and  some  evaluation  of 
several  different  ways  of  accomplishing  step  (c)  above  (con¬ 
fidence  limits  for  the  probability  of  survival  or  excedance 
of  time  h  and  related  topics) .  It  will  discuss  four  differ¬ 
ent  methods  for  attacking  the  estimation  and  confidence 
limits  problem. 

1.  Mathematical  Formulation 

We  shall  assume  that  (x^ ,x^ , . . . ,xn)  are  the  complete 
times  of  repair  (or  down  times) ,  and  that  they  are  indepen¬ 
dent  realizations  of  the  generic  random  variable  X,  where 

2 

Y  »  lnX  is  normally  distributed  with  mean  y  and  variance  a  , 
both  unknown.  This  kind  of  assumption  is  often  made  in  prac¬ 
tice.  This  implies  that  the  probability  that  a  randomly  se¬ 
lected,  future,  down  or  repair  time  exceeds  h  is  given  by  the 
formula 


(1.1 


In  practice,  this  formula  is  not  immediately  applica- 
2 

ble  when  y  and  a  are  unknown  but  if  we  estimate 


c2  -  2  dnxi  -  Inx)2  , 


(1.3) 


vie  can  go  ahead  and  quote  a  point  estimate;  the  latter  de¬ 
pends  on 


lnh  -  jj 

Vs2” 


(1.4) 


If  we  examine  this  quantity  (integration  limit)  it 
will  be  seen  to  be  a  single  realization  of  a  random  variable 
written  as 


(1.5) 


—  2  2  2  2 

where  Y  is  N(y,a  )  and  is  a  X(n_]j  proportional  to  a 

Chi-squared  r.v.,  the  latter  being  independent  of  Y  by  the 
convenient  (log)  normal  assumption.  Now  re-write  0  as 


0  = 


-  (H-y)  ~  (Y-y) 


(1.5) 


/n 


s  2 

_z_ 

n 


If  (H-u)  =  0  then  (-Q/n)  would  be  precisely  distri¬ 
buted  as  a  Student's  t.  On  the  other  hand 


-9/n  = 


(Y-u+  y-H)  /n 


(1.7) 


If  we  write  <$  =  y-H  then 


-e/5  = 

y 


(1.8) 
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has  a  known  density  function,  that  of  the  Non-central  t  which 


is  conveniently  expressed  in  terms  of  the  non-centrality 
parameter 


Classical  methods  exist  for  utilizing  this  to  estab¬ 
lish  tolerance  limits.  In  this  paper  a  different  approach 
is  followed.  We  examine  the  performance  of  several  conve¬ 
nient  approximate  methods  for  assessing  the  uncertainty  in 
the  simple  point  estimate  (1.1),  where  estimate  (1.2)  and 
(1.3)  are  used  for  the  parameter  values.  These  methods  are 
the  Delta  method  -(linearization),  the  Jackknife,  and  the 
Bootstrap,  as  well  as  a  completely  distribution-free  (Ber¬ 
noulli  trials)  method.  Details  now  follow. 


B.  PURPOSE  AND  APPROACH 

1 .  Distribution-Free  Approach 

In  general,  suppose  we  want  to  solve  the  problem  of 
estimating  the  survival  probability  without  any  distribu¬ 
tional  assumption,  other  than  that  observations  are  iid.  The 
simplest  way  is  to  use  the  binomial  approach.  If  (x1,x2,..., 

x  )  indicates  the  iid.  sample  of  down  or  repair  times,  we  can 
n 

estimate  P (x>h) ,  survival  probability,  by  means  of 


P  [X>h]  =  #  ( ——  =  p 
—  n 


(1.9) 


Then  we  can  set  up  a  confidence  interval  for  the  sur¬ 
vival  probability  (1.9)  by  making  use  of  the  fact  that  for 
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rge  n  the  binomial  distribution  can  be  approximated  by  a 
rmal  distribution.  An  approximate  (1-a) *100%  confidence 
terval  for  the  binomial  parameter  p  is  given  by 


p  -  z 


<  P  <  P  + 


(1.10) 


ere  za^2 (1  -  a/2) *100%  point  of  the  tabled  unit  normal. 

2 .  Maximum  Likelihood  Approach 

We  can  assume  along  with  others,  that  repair  time 

ta  comes  as  a  random  sample  from  a  log-normal  population: 

2 

*  lnX  where  X  is  Normal (y, a  ).  This  assumption  will  be 
ucial  in  all  three  methods.  Then  the  maximum  likelihood 
timates  (M.L.E.)  of  the  parameter  are  as  stated  before: 


o  -  y  -  k 

i=l 


(1.11) 


92  -  n^T  £.‘yi  -  ?>2  *  Sy2 


(1.12) 


Strictly  speaking,  a  (1.12)  is  the  M.L.E.  multiplied 

2 

(n/n-1)  and  is  unbiased  for  c  .  Furthermore,  in  repeated 
mples  of  size  n  we  have  that  "exactly"  (assuming  the  model 
correct) : 


ft  *  Y  is  Normal  (u,—) 


(1.13) 


(ii)  d4 


2  .  a2x2(n-l) 

y  13  n=T 


(1.14) 


where  E[d2]  =  a2,  and  Var[d2]  =  2a4/(n-l) 


(iii)  0  and  d *  are  statistically  independent. 
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Thus  for  large  n  both  0  and  a  tend  to  be  close  to 
their  respective  population  values,  guaranteeing  a  good  ap¬ 
proximation  to  the  survival  probability  if  model  (1)  is  cor¬ 


rect.  Now  according  to  the  assumed  model,  the  probability 
of  exceeding  h  hours  is 


P (x>h) 


(1.15) 


The  maximum  likelihood  estimate  of  this  probability 

r*  2  2 

is  obtained  by  replacing  y  by  0,  a  by  d  . 


(1.16) 


Now  find  upper  and  lower  limits  for  the  parameter: 

q  .ISkJLJL  (1 

i.e.  q  and  3  are  functions  of  the  observations  such  that 
£  -  q  -  q  with  prescribed  probability  ( 1— ct)  *100%,  say  95%. 
These  then  translate  into  upper  and  lower  limits  on  the 
probability  of  exceeding  h 


00  ^  00 

(  dz  f  -*Z  dz 

/  e  —  <  P  (x>h)  <  /  e  — 


(1.18) 


If  we  compute  q  and  q  from  a  sample,  then,  under  the 
initial  assumptions,  we  have  the  desired  confidence  limits 
for  the  probability  of  exceeding  h. 

3.  Delta  Method  (DL) 


The  delta  method  is  an  approximate  way  of  finding 
the  distribution  of  q.  It  is  known  that  functions  such  as  q 
are  approximately  normally  distributed  for  "sufficiently 
large"  n  (see  Cramer  [Ref.  2}).  We  estimate  q  by 


lnh  -  fl 
4 - it - 


(1.19) 


and  use  the  "delta  method",  or  method  of  linearization,  or 
small  errors,  to  estimate  the  variance: 


2  2 

Var [§]  s  (||)  VartCl]  Var[d2] 


(1.20) 


There  is  no  covariance  term  because  of  the  (theoretical)  in-  ' 

2 

dependence  of  0  and  d  ,  see  (iii)  above  in  section  2.  This 
formula  yields 


3q  _  1  „ _ r<M_o  _ _ r*2,  _  2a  3q_-(lnh-0) 

9?  “d'  VartB1  Var[4  1  n^l'  ~T72  z  TTrrm  • 


3d*  2 • (d  ) 


(1.21) 


vwni  ,  i  .£  +  i.uqh  -  a>2.  2gj. 

q2  n  4  (a2)3  (n“1) 


Ml  -  i  + 


d  • (n-1) 


r 

1  i  .  1 

— <1  +  =■* 
n  2 


(lnh-Q) 


(1.22) 


(1.23) 
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Assume  3  can  be  taken  to  be  normal  with  mean  q  and 
2 

variance  d  and  quote  these  approximate  confidence  limits: 


■Jdl  ■  4  + 


(1.24) 


2DL  *  «  '  Zl-a/2W 


(1.25) 


This  translates  into  the  desired  (but  approximate) 
confidence  limits  for  the  probability  of  exceeding  h: 


<  P(x>h) 


(1.26) 


Several  approximations  have  been  made  in  the  process 
described  and  the  validity,  for  moderate  n,  of  such  a  rela¬ 
tively  simple  process,  must  be  checked.  Notice  that  the 
exact  distribution  of  q  is  non-central  t  under  the  basic 
model  assumption.  This  approach  replaces  the  n-c.t  by  a 
convenient  normal  approximation. 

4 .  Jackknife  Method  (JK) 

The  jackknife  is  an  alternative  way  of  putting  con¬ 
fidence  limits  on  the  parameter 

lnh  -  u 
a  =  - - 


For  further  discussion  see  Mosteller  and  Tukey  [Ref. 

3]  and  Efron  (Ref.  4].  In  brief,  the  jackknife  method  has 
the  capacity  to  reduce  the  bias  of  estimates  of  such  quantities 
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and,  more  importantly,  to  furnish  confidence  limits  that  be¬ 
have  in  a  satisfactory  manner. 


Jackknife  estimates  and  confidence  limits  are  con¬ 
structed  by  successively  leaving  out  parts  of  the  available 
data  to  construct  pseudovalues.  These  are  then  averaged,  and 
the  stability  of  the  average  assessed  by  use  of  Student's  t 
or  the  Normal  in  order  to  obtain  confidence  limits.  The  pro¬ 
cedure  is  given  below  for  our  case: 

(1)  Form  the  estimate 


qn,yl’y2'y3 


(1.27) 


This  is  the  m.l.e.  using  all  the  data,  just  as  before. 

(2)  Form  the  estimates  q^^j  i *yl'y2' *  *  * 'yi-l ' 
yi' *  *  * ,yn>  i  *  1,  2,  .  .  .  ,  n;  these  are  similar  to  qn, 
but  omit  successively  each  single  observation  yj^  ,y2  *  •  •  •  ,yR; 
at  the  next  stage  each  observation  is  then  restored  and  the 
following  taken  out,  as  i  runs  from  1  to  n  and  thus  there 
are  n  values 

(3)  Compute  the  pseudovalues  as  follows: 

ui  *  nqn  -  (n-l)qtn-1)  ^  i  *  1,  2,  .  .  .  ,  n  (1.28) 

(4)  Compute  the  mean  and  variance  of  the  pseudovalues 


(5)  Approximate  (accuracy  increasing  with  n  in¬ 


creasing)  (1-a) *100%  confidence  limits  for  q  are  given  by: 

s  s 

3jk  1  5  *  ^0/2  (”-1>  7=^/2  (n_1)  5  <l-31> 

where  *3/2  is  the  d  ~a/2)*100  percent  point  of  Student's 

t  (the  standard,  central,  distribution).  Also  we  can  use  za/2 
as  before  as  an  option. 

(6)  This  means  that,  with  approximate  (1-a) *100% 
confidence,  the  probability  of  survival  is  between  the  two 
confidence  limits  that  follow: 


i  P (x>h) 


(1.32) 


This  procedure ,  based  on  the  m . 1 . e . ,  has  been  theo¬ 
retically  validated  for  large  n.  It  competes  with  the  delta 
method,  but  is  somewhat  more  difficult  to  carry  out. 

5 .  Bootstrap  Method  (BT) 

The  bootstrap  method  (see  Efron  (1979)  [Ref.  4J)  is 
similar  to  the  jackknife  method,  but  differs  in  being  a  re¬ 
sampling  procedure.  The  procedure  is  given  as  below  for  our 
case: 


(1)  Calculate 

4  *  <*»'*! '*2 


This  is  the  m.l.e.  using  original  data,  same  as  before. 


(2)  Draw  a  "Bootstrap  sample”,  using  y  ^yj#  •  ••  *Yn  as 
basic  distribution,  value  each  having  probability  1/n 


*1^2 '*3 . *n 


and  calculate  u  =  S^Yi'Y^' '  *  * 'yn* 


lnh  -  y* 


(3)  Independently  repeat  step  (2)  a  large  number  of 
times,  B,  obtaining  "bootstrap  replications"  u^  i  =  1,2,...,B, 
and  calculate 


1  a 
5  -  * 

2  1?  _ 

=  STT  J1  (ui  -  »> 


(1.35) 


(1.36) 


(4)  Approximate  (1-a) *100%  confidence  limits  for  q. 
Here  are  four  different  approaches: 

(a)  Non-Par ame trie  Approach  (BT1) :  Take  the 
order  statistics  of  bootstrap  sample 

u(l)  <  u (2)  <  *  *  *  <  U(B) 

then  let  j  *  £j*bJ,  and  as  confidence  limits  for  q: 

2bT  =  U(J)  "  q  "  U(B-J+1)  =  qBT  U.37) 


(b)  Normal  Approximation  Approach  (BT2) :  If  we 
assume  the  bootstrap  sample  is  approximately  normally  dis¬ 
tributed  then  approximate  confidence  limits  for  q  can  be  set 
down: 


=  u  -  z 


1  -  a/2 *Su  "  q  “  u  +  Z1  -  a/2 'Su  =  qBT 


(1.38) 
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(c)  Bias-Adjusted  Non-Par  suae  trie  Approach  (BT3)  : 
The  bootstrap  estimate  of  bias  is 

BIAS  =  u  -  q 

where  q  is  the  estimate  of  q  from  the  original  data.  Con¬ 


fidence  limits  for  this  case: 


aBT  =  u(j)  -  BIAS  -  -  u(b- 


(B-J+l) 


-  BIAS  i 


(1.39) 


(d)  Bias-Adjusted  Normal  Approximation  Approach 
(«T4) :  In  this  approach  the  confidence  limits  are: 

3bt=«}-bSS)  (4-^)  +z1.a/2-Su=qBT  (1.40) 

Some  simulation  results  for  these  four  approaches 
will  be  presented  in  the  analysis  section.  Identification 
for  these  cases  are  BT1 ,  BT2 ,  BT3,  BT4. 
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II.  SIMULATION  PROCEDURE 


A  simulation  procedure  has  been  used  to  compare  the  three 
methods  for  obtaining  confidence  intervals  for  the  probability 
a  repair  time  exceeds  h  ■  4  hours  in  the  case  in  which  the 
log-normal  assumption  is  met  and  other  cases  in  which  it 
isn't  (the  exponential  and  long-tailed  exponential) .  Speci¬ 
fically,  simulation  has  been  used  to  compute 

(a)  The  actual  coverage  of  the  true  survival  probability 
by  the  confidence  intervals  given  by  the  procedures 
under  study,  when  the  nominal  coverage  is  (l-a)*100%. 

(b)  Measure  of  confidence  interval  size:  the  expected 
width  and  standard  deviation  of  width. 

The  simulation  programs  were  written  in  FORTRAN  IV,  and 
the  simulations  have  been  carried  out  on  the  IBM  3033  at  the 
Naval  Postgraduate  School.  The  Naval  Postgraduate  School 
LLRANDOM  package  was  used,  along  with  the  International  Mathe¬ 
matical  and  Statistical  Library  (IMSL)  random  generator;  1000 
replications  were  used  to  evaluate  each  procedure  in  each 
distributional  situation.  Also  B  *  200  bootstrap  replications 
were  taken  for  each  trial,  four  sample  sizes,  n  »  10,  20,  30, 
40,  and  h  ■  4  hours  and  three  distributions:  Log-normal,  Ex¬ 
ponential,  and  a  Long-tailed  Exponential  were  investigated. 

An  outline  of  the  simulation  procedure  now  follows: 

(A)  Log-normal:  In  this  case  the  basic  variables,  the 
down  or  repair  times,  are  i.i.d.  log-normal  (see  Appendix  B) . 


This  case  was  simulated  for  the  following  "population"  pa¬ 
rameter  values: 


(1) 

U  * 

l. 

„2 

»  1 

(2) 

V  * 

l. 

a2 

*  ln2 

(3) 

y  * 

l. 

„2 

*  31n2 

(4) 

y  = 

l. 

o2 

=  ln2/3 . 

The  procedures  of  the  previous  section  were  used  to  ob¬ 
tain  confidence  intervals  for  the  probability  a  repair  time 
exceeds  h  *  4  hours. 

(B)  Stretched  long-tailed  exponential:  Down  or  repair 
times  come  independently  from  a  stretched  long-tailed  expo¬ 
nential;  see  Appendix  B.  The  simulated  data  was  treated  as 
if  it  was  a  sample  from  a  log-normal  distribution  and  proce¬ 
dures  of  the  previous  section  were  carried  out.  In  this 
simulation  we  used  the  log  transformation  to  tend  to  convert 
the  long-tailed  exponential  observations  towards  normality 
(symmetrize  them) .  The  stretched  long-tailed  exponential 
model  is: 

X  -  AZ (1  +  CZ)  (2.1) 

X  is  stretched  long-tailed  exponential  where  Z  has  the  Ex¬ 
pond)  distribution.  Simulation  was  carried  out  for  A*  3.225 
C  *0.1948.  These  values  were  taken  in  order  to  compare  the 
results  with  Exponential  (A  *  0.22)  case. 

(C)  Exponential  down  or  repair  time:  In  this  situation 
simulations  were  carried  out  for  two  cases.  First  taking  the 


log  of  exponential  down  or  repair  times,  and  treating  these 
as  having  the  normal  distribution;  that  is  treating  the  data 
as  log-normally  distributed.  Second,  taking  the  p  power  of 
data  values,  and  then  treating  the  transformed  values  as  nor¬ 
mally  distributed.  Also  Appendix  C  gives  an  algorithm  for 
estimating  the  p  value  from  the  data.  Simulations  are  also 
carried  out  for  p  *  0.33  (the  classical  Wilson-Hilferty  value) 
and  for  three  X  values. 


.v«. 


III.  ANALYSIS 


The  methods  for  obtaining  95%  confidence  intervals  for 
the  probability  that  a  repair  time  exceeds  h  =  4  hours  de¬ 
scribed  in  chapter  1  were  performed  on  simulated  data  having 
various  distributions;  these  distributions  were  described  in 
chapter  2.  Simulation  results  for  each  method  are  shown  in 
Tables  1  to  11. 

If  we  examine  these  tables  case  by  case,  we  can  find 
these  results: 

(A)  Log-normal  data:  All  three  methods  work  very  well 
for  this  case  except  BT1  and  BT3 ;  they  seem  to  consistently 
have  less  than  nominal  coverage.  The  simple  delta  method 
exhibits  good  coverage,  and  always  has  relatively  small 
average  width,  and  also  low  standard  deviation.  It  seems  to 
work  as  well  as  uX  and  BT4  for  small  sample  sizes  (n  »  10, 

20) .  In  large  sample  sizes  all  methods  agree  in  their  cover¬ 
age  except  for  BT1 .  The  method  BT4  always  appears  to  exhibit 
over-coverage . 

(B)  Stretched  long-tailed  exponential  data:  Table  5 
shows  that  JK  and  BT4  exhibit  over-coverage  when  the  sample 
size  n  *  10.  JK,  DL,  BT4  appear  to  exhibit  correct  coverage 
for  n  »  20,  30,  40,  the  others  don't;  especially  at  sample 
size  n  *  40  they  are  very  poor.  Also  there  is  decreased 
coverage  for  all  methods  when  sample  size  increases;  the 


JK  and  DT4  methods  have  lower  average  width  than  does  DL, 
when  sample  size  increases.  This  is  results  of  the  bias. 

(C)  Exponential  case:  Tables  6,  7,  3  show  that  the  log- 
transformation  may  give  very  poor  results  for  the  exponential 
case,  especially  at  Table  7.  If  we  examine  Table  6  and  Table 
8,  these  tables  show  DL,  JK,  BT4  work  well  for  small  sample 
sizes.  Tables  9,  10  and  11  indicate  that  the  power  transfor¬ 
mation  works  better  than  the  log  transformation.  The  JK, 

BT2 ,  and  BT4  methods  always  have  better  coverage  than  the  DL 
method.  Also  all  methods  agree  in  their  coverage  when  the 
sample  size  increases,  as  was  true  for  the  log-normal  case. 
Generally  JK  and  -BT4  exhibit  acceptable  coverage. 
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bie  1:  Simulation  results  for  log- Normal ( 1 a=1. » 


h  =4.  0 


Samp  le 

A  verags 

St  d . D  9 v 

Siz  = 

Method 

Co  ver  aga 

Width 

Width 

10 

DL 

0 .9420 

0.4353 

0.0409 

JK 

0.9520 

0.4985 

0.  1202 

3?1 

0.8940 

0.4683 

0.0979 

3T2 

0 . 9480 

0.5105 

0.1131 

ST  3 

0 .9150 

0.4857 

0,0940 

BT4 

C .9720 

0.5380 

0. 1205 

20 

DL 

0.9380 

0.3221 

0.0184 

JK 

C . 95  60 

0.3455 

0.0510 

BT  1 

i 

0.9060 

0.3339 

0.0461 

BT2 

0 .9400 

0. 3389 

0.0452 

BT3 

0.9290 

0.3405 

0.0469 

BT4 

0.9620 

0.3495 

0.0491 

30 

DL 

0.9450 

0.2670 

0.0114 

JK 

0 .9520 

0.2795 

0.0302 

BT  1 

0.9240 

0.2748 

0.0306 

BT2 

0.94  30 

0.2737 

0.0296 

BT3 

0  ,93  80 

0.2772 

0.0308 

BT4 

0.9550 

0.2804 

0.0298 

40 

DL 

0.9500 

0.2329 

0. 0084 

JK 

0.9560 

0.2408 

0.0215 

BT  1 

0 .9230 

0.2375 

0.0246 

BT2 

0.9360 

0.2369 

0.0227 

BT3 

0.9460 

0.2385 

0.0234 

BT  4 

0 .96  10 

0. 2406 

0.0226 

case- 
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Table  2:  Simulation  results  for  log-Normal  (ys1 .  ,a=.83)  case 


h=4.  0 


Sample 

Average 

St d. Dev 

Siz  e 

Method 

Coverage 

Width 

Width 

10 

OL 

0.94  23 

0.4278 

0.0468 

JK 

0.9553 

0.4931 

0.1284 

BT 1 

0.8753 

0.4464 

0.1069 

3T2 

0.9393 

0.4952 

0.1151 

BT3 

0.9190 

0.4717 

0. 1002 

BT4 

0.9753 

0.5352 

0.1281 

23 

OL 

0.9370 

0. 3165 

0.0218 

JK 

0.9520 

0.3407 

0.0567 

BT  1 

0.9130 

0.3259 

0.0506 

BT2 

0.9383 

0.3330 

0.0499 

BT3 

0.9283 

0.3325 

0.0506 

BT4 

0.96  00 

0.3452 

0.0538 

33 

DL 

0.94  50 

0.2624 

0.0137 

JK 

0.95  23 

0.2755 

0.0342 

BT  1 

0.9303 

0.2681 

0.0332 

BT2 

0.9430 

0.2704 

0.0316 

BT3 

C .9390 

0.2713 

0.0334 

BT4 

0.9540 

0.2764 

0.0330 

43 

DL 

0 .9520 

0.2289 

0.0101 

JK 

0.95  83 

0.2369 

0.0246 

3T1 

0.9310 

0.2316 

0.0250 

BT2 

0 .9440 

0.2332 

0.0240 

3T3 

0.9470 

0.2335 

0.0253 

BT4 

0.9560 

0.2363 

0.0249 
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J-  y  ~9  i.  *a. 

p  ,  > 

»*  _  . 

Table  3:  Simulation 

results  for 

log-Hormal  (  jj= 

1 . ,a=.48)  case. 

a 

o 

• 

sf 

II 

JS 

Sample 

Average 

Std.Dev 

Size 

Method 

Coverage 

Width 

Width 

i  ■  10 

DL 

0.9430 

0.3823 

0.0700 

1 

JK 

0.9580 

0.4621 

0.1640 

3T1 

0.85  30 

0.  3525 

0.1263 

i/ . 

BT2 

0.9490 

0.4404 

0.1440 

i 

BT3 

0.9190 

0.3944 

0.1218 

BT4 

0  .9810 

0.5180 

0.1615 

V-  23 

DL 

0.9460 

0.2810 

0.0366 

JK 

0.9520 

0.3104 

0.0771 

% 

BT  1 

0.9020 

0.2710 

0.0630 

BT2 

0.9380 

0.2895 

0.0646 

BT3 

0  .9290 

0.2845 

0.0633 

By. 

0.9580 

0.3167 

0.0714 

■ 

*  33 

DL 

0.9460 

0.2325 

0.0241 

JK 

0.9480 

0.248-3  - 

0.0480 

,  - 

BT  1 

0.9220 

0.2275 

0.0417 

i 

BT2 

0.9420 

0.2360 

0.0421 

3T3 

0.9350 

0.2345 

0.0421 

* . ' 

3T4 

0.9510 

0.2500 

0.0446 

43 

DL 

0.9580 

0.2022 

0.0181 

* 

JK 

0.9560 

0.2113 

0.0349 

BT  1 

0.9260 

0.1979 

0.0317 

3T2 

0.9480 

0.2033 

0.0315 

3T3 

0 .94  50 

0.2024 

0.0322 

« 

BT4 

0.9550 

0.2121 

0. 0333 

H  ' 

»  * 

1 

26 

.*• 

»"* 

-*♦  ---*  i»  "  *  *-«-  j  "-* 

-.1 

Table  4:  Simulation  results  for  log-Bormai  (y*  1.  ,a=1. 44) 
case. 

h*4.  0 


Sample 

Average 

Std.Dev 

Size 

Hethod 

Coverage 

Width 

Width 

10 

DL 

0.9420 

0.4440 

0.0324 

JK 

0  .9520 

0.5046 

0.1093 

BT 1 

0.8880 

0.4851 

0.0886 

BT2 

0  .9430 

0.5169 

0.0988 

BT3 

0.9210 

0.5027 

0. 0849 

BT4 

0.9730 

0.5417 

0.1101 

20 

DL 

0.9410 

0.3285 

0.0135 

JK 

0.9550 

0.3503 

0.0423 

BT  1 

0.9180 

0.3464 

0.0410 

BT2 

0.9400 

0.3483 

0.0394 

BT3 

0.9290 

0.3499 

0.0413 

BT4 

0.96  20 

0.3544 

0.0421 

30 

DL 

0.9440 

0.2724 

0.0081 

JK 

0.95  20 

0.2862 

0.0241 

BT1 

0.9280 

0.2825 

0.0271 

BT2 

0.9390 

0.2922 

0.0245 

BT3 

0.9360 

0.2841 

0.0273 

BT4 

0.9560 

0.2849 

0.0254 

40 

DL 

0.95  20 

0.2375 

0.0059 

JK 

0.9600 

0.2451 

0.0165 

BT  1 

0.9390 

0.2435 

0.0209 

3T2 

0.94  80 

0.2434 

0.0186 

3T3 

0.94  90 

0.2444 

0.0210 

BT4 

0.9570 

0.2449 

0.0191 

I 


i 


C-. 


I*  1 
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Table  5:  Simulation  results  for  stretched  Long-tailad 
Exponential. 


h*4.0  4*3.225  C*0.1948 


Sample 

Average 

Std.Dev 

Size 

Method 

Co verage 

Width 

Width 

10 

DL 

0.9530 

0.4341 

0.0350 

JK 

0.9720 

0.4973 

0.1094 

BT1 

0.8550 

0.4404 

0.1052 

BT2 

0 .92  80 

0.4632 

0. 1055 

BT3 

0.9093 

0.4563 

0.1001 

BT4 

0.9853 

0.4893 

0. 1058 

23 

DL 

0.96  50 

0.3195 

0.0170 

JK 

0.95  90 

0.3170 

0.0499 

BT1 

0.8910 

0.3013 

0.0485 

BT2 

0.9170 

0.3015 

0.0460 

BT3 

0.9200 

0.3044 

0.0476 

BT4 

0.95  90 

0.3079 

0.0453 

30 

DL 

0.9570 

0.2642 

0.0110 

JK 

0.9463 

0.2493 

0.0321 

BT1 

0.90  20 

0.2424 

0.0333 

BT2 

0.9153 

0.2415 

0.0313 

BT3 

0.9170 

0.2440 

0.0329 

3T4 

0.94  23 

0.2445 

0.0309 

43 

DL 

0.9430 

0.2301 

0.0087 

JK 

0.9163 

0.2115 

0.0245 

BT1 

0 .8753 

0.2076 

0.0260 

BT2 

0.8940 

0.2063 

0.0244 

BT3 

0.89  50 

0.2085 

0.0258 

BT4 

0.9243 

0.2087 

0.0241 

23 


s%ss  using  Log 
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Table  6:  Siaulation  results  for  EXP(X) 
transformation. 


Sample 

h=4.0 

X-0.22 

Size 

Method 

Coverage 

10 

DL 

0.9480 

JK 

0.9  290 

ET1 

0.3450 

BT2 

0.9090 

BT3 

0.9010 

B?4 

0.9740 

20 

DL 

0.9520 

JK 

0.9060 

BT1 

0.8680 

BT2 

0.8840 

BT3 

0.8930 

BT4 

0.9310 

30 

DL 

0.9380 

JK 

0.8790 

BT1 

0.8530 

BT2 

0.8680 

BT3 

0.8810 

e* 

CO 

0.8990 

40 

DL 

0.9010 

JK 

0.3  280 

BT1 

0.8  110 

BT2 

0.8  160 

BT3 

0.8250 

BT4 

0.8440 

P  (X>4. 0)  *0.4148 


Average 

3rd. Dev 

Point 

Width 

Width 

Estimate 

0.4432 

0.0265 

0.  3731 

0.4476 

0. 1000 

0.  3689 

0.4650 

0.0989 

0.  3707 

0.4808 

0.1043 

0.3707 

0.4772 

0.0974 

0.3707 

0.4975 

0.1036 

0.3763 

0.3261 

0. 0127 

0. 3685 

0.3025 

0. 0473 

0.  3665 

0.3133 

0.0489 

0.3670 

0.3105 

0.0457 

0.3670 

0.3151 

0.0483 

0.  3670 

0.3136 

0. 0445 

0.3702 

0.2697 

0.0031 

0.  3661 

0.2416 

0.0310 

0.  3651 

0.2497 

0.0338 

0.  3649 

0.2467 

0.0311 

0.  3649 

0.2504 

0.0335 

0.  3649 

0.2481 

0.0306 

0.3674 

0.2349 

0. 0065 

0.  3650 

0.2068 

0.0237 

0.  3644 

0.2124 

0. 0265 

0. 3641 

0.2105 

0.0245 

0.3641 

0.2128 

0. 0263 

0. 3641 

0.2112 

0. 0241 

0. 3660 
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Table  7:  Siaulation  results  for 

EIP( X)  za 

se  using 

Log 

+ 

w 

ransforaati 

on. 

h*4.0 

X*0 .13 

P  (X>4. 0) 

*0.5945 

Saaple 

average 

Std. Dev 

Point 

Size 

Hethod 

Coverage 

Width 

Width 

Sstiaate 

10 

DL 

0.9100 

0.4475 

0.0298 

0.  5492 

JK 

0.9010 

0.5282 

0.1199 

0.5218 

BT1 

0.8690 

0.5064 

0.0922 

0.5718 

BT2 

0.9110 

0.5490 

0.1102 

0. 5718 

BT3 

0.9020 

0.5280 

0.0932 

0.5718 

BT4 

0.9410 

0.5728 

0.1213 

0.  5245 

20 

DI 

0.8810 

0.3329 

0.0103 

0.5385 

JK 

0.8660 

0.3774 

0. 0665 

0.5250 

BT1 

0.8  970 

0.3705 

0.0462 

0.5500 

BT2 

0.8820 

0.3740 

0.0521 

0.5500 

BT3 

0.9010 

0.3750 

0.  0502 

0.5500 

BT4 

0.8840 

0.3796 

0.0585 

0.5267 

30 

DL 

0.8470 

0.2763 

0.0057 

0. 5343 

JK 

0.8390 

0.3053 

0.0409 

0.5  256 

BT1 

0.8  840 

0.3031 

0.  0319 

0. 5421 

BT2 

0.8600 

0.3018 

0.0322 

0.  5421 

BT3 

0.8890 

0.3048 

0.0336 

0. 5421 

BT4 

0.8460 

0.3043 

0.0354 

0.5264 

«0 

DL 

0.7900 

0.2411 

0. 0035 

0.5315 

JK 

0.7810 

0.2615 

0.0293 

0.  5251 

BT1 

0.8  330 

0.2594 

0.0266 

0.5375 

BT2 

0.8  120 

0.2594 

0.0256 

0.5375 

BT3 

0.8  370 

0.2604 

0.0276 

0. 5375 

BT4 

0.7860 

0.2607 

0. 0273 

0.  5255 

30 


Table  8:  Simulation  results  for 
transformation. 


EXP (A)  rase  using  Log 


h*4.0 

A  *0 .26 

P  (X>4. 0) 

=0.3535 

Sample 

Average 

Std.Dev 

Print 

Size 

Method 

Coverage 

Width 

Width 

Estimate 

ID 

PL 

0.9630 

0.4322 

0. 0343 

0.3197 

JK 

0.9530 

0.4188 

0. 0994 

0.3225 

BT1 

0.8460 

0.4276 

0.1072 

0.  3097 

BT2 

0.9260 

0.4478 

0.  1062 

0.3097 

BT3 

0.9060 

0.4427 

0.1021 

0.  3097 

BT4 

0.9880 

0.4726 

0.1052 

0.3311 

20 

DL 

0.9660 

0.3173 

0.0171 

0. 3178 

JK 

0.9380 

0.2789 

0.0458 

0.3192 

BT1 

0.8790 

0.2959 

0.0489 

0.  3124 

BT2 

0.8980 

0.2857 

0.0459 

0.3124 

BT3 

0.9070 

0.2841 

0.0473 

0.  3124 

8T4 

0.9  500 

0.2917 

0.0447 

0.3232 

30 

DL 

0.9550 

0.2621 

0.0111 

0.  3161 

JK 

0.9160 

0.2228 

0. 0298 

0.3174 

B71 

0.8700 

0.2282 

0. 0332 

0. 3122 

BT2 

0.8890 

0.2271 

0. 0309 

0.3122 

BT3 

0.8  960 

0.2298 

0.0327 

0.  3122 

BT4 

0.9  350 

0.2301 

0. 0302 

0.3199 

40 

DL 

0.9370 

0.2282 

0. 0088 

0.3155 

JK 

0.8840 

0.1910 

0.0229 

0. 3167 

BT1 

0.8  390 

0.1947 

0.0258 

0. 3127 

3T2 

0.8510 

0.1938 

0.0241 

0. 3127 

BT3 

0.8650 

0.1956 

0.0255 

0. 3127 

BT4 

0.8  940 

0.1956 

0.0236 

0. 3185 
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Table  9:  Simulation  results  for  EXP(X)  rase  using  X 
transformation. 


h*4. 

0  X*0.  22 

p*0. 33 

Sample 

Average 

Std.  Dev 

Siz  3 

Method 

Co verage 

Width 

Width 

10 

DL 

0.9300 

0.4460 

0.0302 

JK 

0.9410 

0.5038 

0.0944 

BT 1 

0.86  90 

0.4941 

0.0857 

BT2 

0.9330 

0.5234 

0.0929 

3T3 

0.9050 

0.5101 

0.0829 

BT4 

0.9760 

0.5459 

0.0990 

20 

DL 

0.9390 

0.3307 

0.0118 

JK 

0.9610 

0.3555 

0.0363 

BT1 

0.9150 

0.3515 

0.0356 

BT2 

0.9450 

0.3535 

0.0335 

BT3 

0.9350 

0.3545 

0.0359 

BT4 

0.96  30 

0.3588 

0.0360 

30 

DL 

0.9410 

0.2742 

0.0069 

JK 

0.9570 

0.2883 

0.0200 

BT1 

0.9210 

0.2875 

0.0250 

BT2 

0.9420 

0.2874 

0.0220 

BT3 

0.9310 

0.2887 

0.0254 

3T4 

0.9570 

0.2897 

0.0232 

40 

DL 

0.9270 

0.2390 

0.0055 

JK 

0.9440 

0.2491 

0.0133 

3T1 

0.9200 

0.2480 

0.0193 

3T2 

0.9370 

0.2479 

0.0163 

BT3 

0.9270 

0.2487 

0.0194 

BT4 

0.9440 

0.2492 

0.0167 

Table  13:  Simulation  results  for  EXP(x)  case  using  Xp 
transfor mation. 


h=4.0  X  *3.  13  p=0 .33 


Sample 

Average 

Std.Dev 

Size 

Sezhod 

Coverage 

Width 

Width 

13 

DL 

0.9353 

0.4443 

0.0291 

JK 

0.9533 

0.4911 

0.0912 

BT 1 

0.86  73 

0.4804 

0.0896 

BT2 

0.9493 

0.5133 

0.0933 

BT3 

0.9370 

0.4963 

0.0843 

3T4 

0.9763 

0.5377 

0. 1005 

23 

DL 

0.94  80 

0.3282 

0.0127 

JK 

0.96  70 

0.3411 

0.0351 

BT  1 

0.9153 

0.3409 

0.0403 

BT2 

0.9433 

0.3427 

0.0372 

BT3 

0.9353 

0.3439 

0.0403 

BT4 

0.9703 

0.3479 

0.0385 

33 

DL 

0.9513 

0.2723 

0. 0084 

JK 

0.9613 

0.2762 

0.0195 

BT1 

0.9183 

0.2775 

0.0274 

BT2 

0.9450 

0.2773 

0.0235 

3T3 

0.9360 

0.2783 

0.0273 

BT4 

0.9710 

0.2792 

0.0237 

43 

DL 

0 .94  20 

0.2372 

0.0060 

JK 

0.9520 

0.2383 

0.0140 

BT1 

0.9133 

3.2395 

0.0213 

BT2 

0.9373 

0.2386 

0.0173 

BT3 

0.9223 

0.2392 

0.0212 

BT4 

0.9530 

0.2399 

0.0179 
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Table  11;  Simulation  results  for  EXP(X)  case  using  X 
transformation. 


h=4.0  X=0.26  p=0.  33 


Sample 

Average 

Std. Dev 

Size 

Hethod 

Cc  ver age 

Width 

Width 

ID 

DL 

0.9270 

0.4345 

0.0420 

JK 

0.9390 

0.4970 

0.1106 

BT 1 

0.8660 

0.4677 

0.1027 

BT2 

0.9310 

0.5082 

0.1072 

BT3 

0.9100 

0.4892 

0.0974 

BT4 

0.9760 

0.5415 

0.1127 

2D 

DL 

0.9380 

0.3224 

0.0136 

JK 

0 .9570 

0.3505 

0.0482 

3T1 

0.9D6D 

0.3381 

0.0436 

3T2 

0.9423 

0.3439 

0.0423 

BT3 

0.92  90 

0.3434 

0.0435 

BT4 

0.96  10 

0.3534 

0.0453 

3D 

DL 

0.9380 

0.2672 

0.0117 

JK 

0 .9503 

0.2345 

0.0285 

BT1 

0.9220 

0.2781 

0.0294 

BT2 

0.9430 

0.2805 

0.0277 

BT3 

0.9290 

0.2805 

0.0297 

3T4 

0.9530 

0.2851 

0.0291 

4D 

DL 

0.9240 

0.2323 

0. 0094 

JK 

0.9390 

0.2454 

0.0198 

BT1 

0.9150 

0.2412 

0.0225 

BT2 

0  .9290 

0.2422 

0. 0207 

BT3 

0.9210 

0.2427 

0.0226 

3T4 

0.9410 

0.2451 

0.0212 

IV.  EXAMPLE i  APPLICATION  TO  OPERATIONAL  DATA 


In  this  chapter  four  methods  were  applied  to  a  real  data 
set.  The  methods  are  Binomial  (BN) ,  Delta  (DL) ,  Jackknife 
(JK) ,  Bootstrap  (see  section  I-b) ;  specifically,  the  data 
refer  to  recovery  times  from  loss  of  offsite  power  at  nuclear 
plants.  The  problem  was  to  estimate  survival  probabilities 
for  h  =  1.5,  2.0,  2.5,  3.0,  3.5,  4.0  (hours).  Data  points 
(n  *  42)  are  shown  in  appendix  D.  We  initially  applied 
several  statistical  goodness-of-fit  tests  to  inquire  into 
the  evidence  for  the  adequacy  of  the  Log-Normal  distribution 
as  a  model  for  these  data.  The  results  of  these  goodness  of 
fit  tests  are  as  follows: 

(1)  Chi-square  test:  See  Arnold,  D.  [Ref.  5];  this  ac¬ 
cepts  the  log-normal  model  at  the  significance  level 
a  =*  0.05. 

(2)  Kolmogorov- Smirnov  test:  See  Arnold,  D.  [Ref.  5]; 
this  test  rejects  the  log-normal  model  with  the  tabu¬ 
lated  value  C  *  0.1367  and  test  statistic  D  =  0.21  for 
a  -  0.05. 

(3)  Wilk-Shapiro  test:  See  Hahn,  G.  J.  [Ref.  6];  this 
test  accepts  log-normal  model  for  a  *  0.05. 

We  applied  four  estimation  methods  to  these  data,  utiliz¬ 
ing  the  log-normal  assumption.  The  results  are  shown  in 


Table  12. 


Table  12:  Recovery  Time  Example  Results  (h  in  hours) 


H 


1.5 


2.0 


2.5 


3.0 


Upper  Conf. 

Po  int 

Lower  Conf. 

Method 

Limits 

Esti nation 

Limits 

Width 

BN 

0.476 

0.  333 

0.  191 

0.285 

DL 

0.445 

0.  324 

0.219 

0.226 

JK 

0.432 

0.  324 

0.230 

0.202 

BT1 

0.435 

0.  318 

0.215 

0.220 

BT2 

0.438 

0.  318 

0.214 

0.224 

BT3 

0.441 

0.  318 

0.215 

0.226 

BT4 

0.451 

0.  330 

0.  224 

0.227 

BN 

0.422 

0.  286 

3.  149 

0.  273 

DL 

0.400 

0.  280 

0.182 

0.218 

JK 

0.334 

3.  282 

0.  195 

0.  189 

BT1 

0.379 

3.  278 

0.  198 

0.  180 

BT2 

0.374 

0.  278 

0.  197 

0.  177 

BT3 

0.381 

3.  278 

0.198 

0.  182 

BT4 

0.378 

3.  282 

0.200 

0.  178 

3N 

0.422 

0.  286 

0.  149 

0.273 

DL 

0.366 

3.  249 

3.  155 

0.  211 

JK 

0.348 

3.  250 

0.  169 

0.  179 

BT1 

0.330 

0.  250 

0.  175 

0.  155 

BT2 

0.345 

0.  250 

0.171 

0.  174 

BT3 

0.329 

0.  250 

0.  175 

0.  153 

BT4 

0.342 

3.  247 

0.  169 

0.  173 

BN 

0.395 

3.  26  2 

3.  129 

0.  266 

DL 

0.339 

0.  224 

0.  135 

0.204 

JK 

0.320 

3.  226 

0.  150 

0.  170 

BT1 

0.309 

3.  223 

0.  135 

0.  174 

BT2 

0.324 

0.  223 

3.  143 

0.  182 

BT3 

0.311 

3.  223 

0.135 

0.  176 

BT4 

0.327 

3.  225 

0.  144 

0.  183 

f 


t 

J 

•  i 


4.0 


1 

I 

i 

\ 

t 

> 

i 

\ 

\ 

\ 


BN 

0.395 

0.  262 

0.  129 

0.  266 

DL 

0.3  18 

0.  204 

0.  120 

0.  193 

JK 

0.298 

0.  207 

0.  1  35 

0.  1  63 

BT  1 

0.287 

0.  204 

0.  140 

0.  147 

BT2 

0.239 

0.  204 

0.136 

0.  152 

BT3 

0.288 

0.  204 

0.  140 

0.  147 

ST  4 

0,290 

0.  205 

0.137 

0.153 

3N 

0.367 

0.  238 

0.  109 

0.  258 

DL 

0.299 

0.  188 

0.  107 

0.  1  92 

JK 

0.279 

0.  191 

0.  122 

0.  157 

3T1 

0.263 

3.  185 

0.113 

0.  150 

BT2 

0.273 

0.  185 

0.  1  17 

0.  156 

BT  3 

0.267 

0.  185 

0.  113 

0.  1  54 

BT4 

i 

0.292 

0.  192 

0.  122 

0.  160 

i 


: 
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V.  CONCLUSIONS 


The  Delta,  Jackknife  and  Bootstrap  methods  applied  to  the 
log-normal  model  work  well  when  down  or  repair  times  are 
truly  log-normal.  Especially  notice  that  DL,  JK,  and  BT4 
seem  to  work  much  better  than  BT1 ,  BT2 ,  BT3 .  Recall  that 
these  procedures  do  not  appear  sensitive  to  the  population 
variance;  see  section  III.  It  is  comparatively  easy  to  use 
JK  and  3T4  when  sample  size  is  small  (n  =  10,  20).  The  delta 
method  is  always  convenient,  but  especially  when  the  sample 
is  large  (n  =  40  or  more)  because  it  is  a  very  simple  pro¬ 
cedure  to  apply,  requiring  much  less  computation  than  the 
others.  As  Table  12  shows,  the  Binomial  method  gives  some 
idea  of  the  survival  probability  for  practical  purposes. 

Note  that  Binomial  confidence  limits  are  much  wider  than 
those  that  assume  the  log-normal  model. 

Use  of  the  log  transform  on  exponential  data  produces 
biased  estimates  of  survival.  Use  of  the  power  transforma¬ 
tion  with  (p  *  1/3)  always  gives  a  better  coverage  of  the 
survival  probability  when  data  are  exponential.  One  proce¬ 
dure  was  described  in  Appendix  C  for  estimating  the  p  value 
from  data.  Table  13  gives  simulation  results  for  the  expo¬ 
nential  case.  As  our  results  show,,  this  procedure  is  not 
estimating  p  value  correctly.  If  this  procedure  were  to  work 
correctly  (if  it  could  be  calibrated)  then  we  could  use  the 
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transformation  (for  converting  data  towards  the  normal 
>rm)  without  making  any  assumption  (e.g.  this  data  coming 
:om  exponential  or  gamma  or  log-normal,  etc.).  Then  after 
lis  is  done,  methods  DL,  JK,  BT4  might  produce  considerably 
stter  confidence  limits  for  the  actual  survival  probability. 
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APPENDIX  A 


COMPUTER  PROGRAMS 

Simulation  programs  consist  of  two  main  programs  for 
three  methods  (DL,  JK,  BT) .  These  main  programs  compute 
survival  probability  confidence  limits,  and  scores  the  cover¬ 
age  for  each  replication.  Then,  after  1000  replications  the 
program  computes  the  statistics  of  these  parameters  and 
prints  out  the  results. 

There  is  another  program,  called  SURVP.  This  program 
computes  point  estimates,  confidence  limits,  and  widths  of 
confidence  limits  on  survival  probability,  using  the  BN,  DL, 
JK,  BT  procedures  on  a  given  data  set,  under  the  log-normal 
model  assumption. 

Variables  List: 

R  *  Down  or  repair  times. 

R1  =  Log  of  down  or  repair  times. 

RBAR  *  Mean  of  down  or  repair  times. 

RSD  *  Standard  deviation  of  down  or  repair  times. 

GHAT  =  Point  estimation  of  q  parameter  (see  1.17)  for  delta 
method 

GJK  =  Point  estimation  of  q  parameter  for  jackknife  method. 

GBOOT  =  Point  estimation  of  q  parameter  for  bootstrap 
method. 

VARG  *  Variance  of  point  estimation  for  delta  method. 


SE  *  Standard  error  of  point  estimation  for  jackknife 
method . 

PHAT  *  Point  estimation  of  survival  probability. 

BUP  a  Upper  confidence  limit  estimation  of  q  parameter. 

BLOW  *  Lower  confidence  limit  estimation  of  q  parameter. 

CUP  »  Upper  confidence  limit  estimation  of  survival 
probability. 

CLOW  =  Lower  confidence  limit  estimation  of  survival 
probability. 

AINT  =  Width  of  estimated  confidence  limits  of  survival 
probability. 

G  =  Pseudovalues. 

F  *  Bootstrap  replications. 

N  »  Number  of  data  points. 

N1  *  Number  of  replications. 

N2  *  Number  of  bootstrap  replications. 
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APPENDIX  B 


LOG-NORMAL  AND  STRETCHED  LONG-TAILED 
EXPONENTIAL  DISTRIBUTION 


(1)  Let  x  be  log-normal  random  variable  which  ln(x) 
2 

i,a  )  k  moment  of  x  as  follows: 


E[xk]  =  exp(ku  +  l/2*k2a2) 

s  [Ref.  7]  so 

E [x]  =  exp(y  +  l/2*a2) 
E[x2]  «  exp(2y  +  2*a2) 

Var [x]  =  (E [x] ) 2  (e°2  -  1) 
"coefficient  of  variation"  as  follows: 

vsrtxl  .  e3’  .  j 

/  -c*  r  ..  i  \  ^ 


(2)  Let  W  be  stretched  long-tailed  exponential  variable 
Lch  W  *  A* z ( 1  +  C*z)  where  z  is  unit  exponential  and  A  and 
ire  constant.  If  we  write  CDF  for  this  distribution  as: 

P(Wiw)  =  P  [  A*  z  ( 1  +  C  *z)  -  w]  =  P  [  Z— z  (w)  ] 


=  A* z ( 1  +  C*z) .  If  we  solve  this  equation  for  z  we  can  get 
/)  as 


z  (w) 


-A 


t  /a2 

2  AC 


+  4ACw 
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E  [w]  =  A  (1  +  2C) 

E [w2 ]  =  A2  [2  +  12C  +  24C] 

Var [w]  =  A2[l  +  8C  +  20C] 

(cv)2  =  Var  M  =  1  +  8C  +  20C2 
(E [w] ) 2  1  +  4C  +  4C2 

2 

If  we  look  at  w  as  a  log-normal  (u,a  )  variable  then 


1  +  8C  +  200  .a2 

-  -  e  -  i 

1  +  4C  +  4C 


We  can  get  C  value  from  this  equation,  then 


:[w]  =  A  ( 1  +  2C)  =  eu  + 


gives  A  value. 
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APPENDIX  C 


POWER  TRANSFORMATION 


The  problem  is  in  the  x*5  transformation  (toward  the  nor¬ 
mal  form)  finding  the  p  value  for  given  data.  One  method 
for  finding  the  p  value  is  as  follows  [Ref.  8] : 
x1,x2,...,xn  data  points  and  M  »  Median  of  this  data. 

(1)  Take  order  statistic  of  given  data 

.  .  <  x 


X(l)  <  X  (2)  < 


(2)  Then  compute  values  as 


(n) 


1  -  Pj 


(x ln-i+1)  -M)-(M-Xm)  ' 
l_(x(n-j.l)  ‘M)2+(M'x(jl)2J 


2  *M  j  *  1,2, . . . ,n/2 


where  M  is  the  median  of  the  data. 


(3)  Take  the  median  or  mean  of  q ^ : 


q  *  median (q^)  or  q 


1  n/2 
TnTTT  j 


(4)  Then  we  can  get  p  value  as  p  ■  1  -  q. 

Table  13  gives  simulation  results  for  this  algorithm  for 
the  exponential  case.  For  the  exponential  case  the  best  p 
value  is  p  ■  1/3.  Simulation  results  do  not  give  this  value 
so  this  algorithm  is  not  working  correctly. 


Table  13:  Simulation  results  for  ?  algorithm  using 
Exponential  data. 


X»0.22 


Sample 

Average 

Variance 

Skewness 

Size 

of  P 

of  P 

of  P 

10 

0.74  34 

0.9174 

0.7535 

20 

0.6039 

0.3507 

1 .3270 

30 

0.5336 

0.24  55 

1.3819 

40 

0.5121 

0.1872 

1.8234 

50 

0.4903 

0.1397 

1.4981 

60 

0.4859 

0.1122 

1.4225 

70 

0.4658 

0.0912 

1.5475 

80 

0.4571 

0.0814 

1.6633 

90 

0.44  98 

0.0710 

1.5987 

100 

0.4470 

0.0669 

1.8234 

110 

0.4403 

0.05  80 

1.9510 

120 

0.4323 

0.0515 

1.4446 

130 

0.4298 

0.04  34 

1.3154 

140 

0.4263 

0.0452 

1  .4466 

150 

0 . 4 1 96 

0.03  94 

1.5704 

i 
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APPENDIX  D 


MEAN  SQUARE  ERROR  OF  SURVIVAL  PROBABILITY 

Mean-square  errors  were  calculated  for  exponential  case 
using  the  log  transformation  for  three  values  and  h  =  4.0. 
The  procedure  is  as  follows: 

(1)  Generate  exponential  sample  x^,x2,...,xn  , 
n  »  10,20, .. . ,250. 

(2)  Find  actual  survival  probability  as: 

P (x>h)  -  e“Xh 

(3)  Estimate  survival  probability  (incorrectly)  as: 

P(x>h)  * 

lnh 

5" 

(4)  Calculate  (P(h)  -  P(h))2 

(5)  Repeat  this  procedure  1000  times. 

(6)  Calculate  MS^  as: 

/V  ,  1000  - 

MSE  =  l  [P  (h)  -  Pi(h)  ]2 

i»l 

Simulation  results  for  mean  square  error  of  survival 
probability  were  shown  in  Tables  14,  15  and  16. 


Tibia  14:  Mean  Square  Error  of  Survival  probabilit 
for  EXP (X)  case  using  Log  transformation 


P 


Sample 

h*4.0 

Bean  Square 

Size 

Error  (BSE) 

10 

0.0166 

20 

0.0085 

30 

0.0063 

40 

0.0056 

50 

0.0049 

60 

0.0045 

70 

0.0041 

80 

0. 0040 

90 

0.0038 

100 

0.0038 

110 

0.0037 

120 

0.0035 

130 

0.0034 

140 

0.0033 

150 

0.0033 

160 

0.0033 

170 

0.0032 

180 

0.0031 

190 

0.  0032 

200 

0.0032 

210 

0.0031 

220 

0. 0031 

230 

0.0031 

240 

0.0031 

250 

0.0031 

Square  Root 
Of  BSE 

0.1288 
0.0920 
0.0795 
0.0746 
0.0699 
0.0673 
0.0642 
0.0635 
0.0614 
0.0612 
0.0604 
0.0594 
0.0585 
0.0578 
0.0575 
0.0  572 
0.0564 
0.0560 
0.0562 
0.0  563 
0.0561 
0.0557 
0.0558 
0.0557 
0.0553 
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Table  15:  Kean  Square  Error  of  Survival  probability 
for  EXP (X)  case  using  log  transformation. 


h*4.0  A-0.13 


Samp le 

Kean  Square 

Square  Root 

Size 

Error  (BSE) 

of  KSE 

10 

0.  0234 

0.1531 

20 

0. 0128 

0. 1 132 

30 

0. 0098 

0.0991 

40 

0.0089 

0.0942 

50 

0.0081 

0.0897 

60 

0.0076 

0.  0869 

70 

0.0070 

0.0834 

80 

0.0068 

0.0824 

90 

0.0065 

0.0803 

100 

0.0064 

0.0802 

110 

0.0063 

0.0792 

120 

0.0061 

0.0781 

130 

0.0059 

0.0770 

140 

0.0059 

0.0766 

150 

0.0058 

0.0762 

160 

0.0058 

0.0758 

170 

0.0057 

0.0753 

180 

0.0055 

0.0745 

190 

0.0056 

0.0746 

200 

0.0056 

0.0747 

210 

0.  0056 

0.0747 

220 

0.0055 

0.0744 

230 

0.0056 

0.0745 

240 

0.0055 

0.0743 

250 

0.0054 

0.0737 

Table  16:  Mean  Square  Error  of  Survival  probability 
for  EXP (X)  case  using  Log  transformation. 


h=4 . 0  A  =0 

.25 

Sample 

Mean  Square 

Square  B03 

Size 

Error  (MSE) 

of  MSE 

10 

0.0135 

0.1 160 

20 

0.0066 

0.0812 

30 

0.0047 

0.0688 

40 

0.  0041 

0.0637 

50 

0.  0034 

0.0585 

60 

0.0031 

0.0559 

70 

0.0028 

0.0525 

80 

0.0027 

0.0519 

90 

0.0025 

0.0495 

100 

0.0024 

0.0492 

110 

0.0023 

0.0  483 

120 

0.0022 

0.0472 

130 

0.0022 

0.0465 

140 

0.0021 

0.0456 

150 

0. 0021 

0.0453 

160 

0.0020 

0.0448 

170 

0.0019 

0.0439 

180 

0.0019 

0.0435 

190 

0.0019 

0.0437 

200 

0.0019 

0.0438 

210 

0.0019 

0.0435 

220 

0.0019 

0.0431 

230 

0.0019 

0.0431 

240 

0.0018 

0.0429 

250 

0.0018 

0.0425 

60 


APPENDIX  E 


DATA  POINTS  FOR  EXAMPLE 


These  are 

plants: 

recovery 

times  (hours) 

from  LOSP  at 

nuclear 

24.6160 

25.6660 

11.0830 

0.0038 

0.3333 

0.6166 

1.5000 

1.1833 

0.0333 

0.0500 

0.2666 

5.8000 

4.9830 

1.8330 

0.5000 

6.4660 

0.2833 

1.0000 

0.9000 

0.1666 

0.6666 

0.4333 

0.1333 

0.0166 

5.5833 

0.4833 

0.0028 

0.9333 

0.2333 

4.9833 

0.1500 

2.6660 

4.7500 

0.1333 

1.0166 

0.0666 

0.0166 

8.9000 

3.5000 

0.9166 

0.3333 


0.0041 
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